Membrainy: a ‘smart’, unified membrane analysis tool

BACKGROUND: The study of biological membranes using Molecular Dynamics has become an increasingly popular means by which to investigate the interactions of proteins, peptides and potentials with lipid bilayers. These interactions often result in changes to the properties of the lipids which can modify the behaviour of the membrane. Membrainy is a unified membrane analysis tool that contains a broad spectrum of analytical techniques to enable: measurement of acyl chain order parameters; presentation of 2D surface and thickness maps; determination of lateral and axial headgroup orientations; measurement of bilayer and leaflet thickness; analysis of the annular shell surrounding membrane-embedded objects; quantification of gel percentage; time evolution of the transmembrane voltage; area per lipid calculations; and quantification of lipid mixing/demixing entropy. RESULTS: Each analytical component within Membrainy has been tested on a variety of lipid bilayer systems and was found to be either comparable to or an improvement upon existing software. For the analytical techniques that have no direct comparable software, our results were confirmed with experimental data. CONCLUSIONS: Membrainy is a user-friendly, intelligent membrane analysis tool that automatically interprets a variety of input formats and force fields, is compatible with both single and double bilayers, and capable of handling asymmetric bilayers and lipid flip-flopping. Membrainy has been designed for ease of use, requiring no installation or configuration and minimal user-input to operate

Polarizable Water Model for the Coarse-Grained MARTINI Force Field

Coarse-grained (CG) simulations have become an essential tool to study a large variety of biomolecular processes, exploring temporal and spatial scales inaccessible to traditional models of atomistic resolution. One of the major simplifications of CG models is the representation of the solvent, which is either implicit or modeled explicitly as a van der Waals particle. The effect of polarization, and thus a proper screening of interactions depending on the local environment, is absent. Given the important role of water as a ubiquitous solvent in biological systems, its treatment is crucial to the properties derived from simulation studies. Here, we parameterize a polarizable coarse-grained water model to be used in combination with the CG MARTINI force field. Using a three-bead model to represent four water molecules, we show that the orientational polarizability of real water can be effectively accounted for. This has the consequence that the dielectric screening of bulk water is reproduced. At the same time, we parameterized our new water model such that bulk water density and oil/water partitioning data remain at the same level of accuracy as for the standard MARTINI force field. We apply the new model to two cases for which current CG force fields are inadequate. First, we address the transport of ions across a lipid membrane. The computed potential of mean force shows that the ions now naturally feel the change in dielectric medium when moving from the high dielectric aqueous phase toward the low dielectric membrane interior. In the second application we consider the electroporation process of both an oil slab and a lipid bilayer. The electrostatic field drives the formation of water filled pores in both cases, following a similar mechanism as seen with atomistically detailed models

Perforin Rapidly Induces Plasma Membrane Phospholipid Flip-Flop

The cytotoxic cell granule secretory pathway is essential for host defense. This pathway is fundamentally a form of intracellular protein delivery where granule proteases (granzymes) from cytotoxic lymphocytes are thought to diffuse through barrel stave pores generated in the plasma membrane of the target cell by the pore forming protein perforin (PFN) and mediate apoptotic as well as additional biological effects. While recent electron microscopy and structural analyses indicate that recombinant PFN oligomerizes to form pores containing 20 monomers (20 nm) when applied to liposomal membranes, these pores are not observed by propidium iodide uptake in target cells. Instead, concentrations of human PFN that encourage granzyme-mediated apoptosis are associated with pore structures that unexpectedly favor phosphatidylserine flip-flop measured by Annexin-V and Lactadherin. Efforts that reduce PFN mediated Ca influx in targets did not reduce Annexin-V reactivity. Antigen specific mouse CD8 cells initiate a similar rapid flip-flop in target cells. A lipid that augments plasma membrane curvature as well as cholesterol depletion in target cells enhance flip-flop. Annexin-V staining highly correlated with apoptosis after Granzyme B (GzmB) treatment. We propose the structures that PFN oligomers form in the membrane bilayer may include arcs previously observed by electron microscopy and that these unusual structures represent an incomplete mixture of plasma membrane lipid and PFN oligomers that may act as a flexible gateway for GzmB to translocate across the bilayer to the cytosolic leaflet of target cells

Purification and Activity Testing of the Full-Length YycFGHI Proteins of Staphylococcus aureus

Background: The YycFG two-component regulatory system (TCS) of Staphylococcus aureus represents the only essential TCS that is almost ubiquitously distributed in Gram-positive bacteria with a low G+C-content. YycG (WalK/VicK) is a sensor histidine-kinase and YycF (WalR/VicR) is the cognate response regulator. Both proteins play an important role in the biosynthesis of the cell envelope and mutations in these proteins have been involved in development of vancomycin and daptomycin resistance. Methodology/Principal Findings: Here we present high yield expression and purification of the full-length YycG and YycF proteins as well as of the auxiliary proteins YycH and YycI of Staphylococcus aureus. Activity tests of the YycG kinase and a mutated version, that harbours an Y306N exchange in its cytoplasmic PAS domain, in a detergent-micelle-model and a phosholipid-liposome-model showed kinase activity (autophosphorylation and phosphoryl group transfer to YycF) only in the presence of elevated concentrations of alkali salts. A direct comparison of the activity of the kinases in the liposomemodel indicated a higher activity of the mutated YycG kinase. Further experiments indicated that YycG responds to fluidity changes in its microenvironment. Conclusions/Significance: The combination of high yield expression, purification and activity testing of membrane and membrane-associated proteins provides an excellent experimental basis for further protein-protein interaction studies an

Determining Peptide Partitioning Properties via Computer Simulation

The transfer of polypeptide segments into lipid bilayers to form transmembrane helices represents the crucial first step in cellular membrane protein folding and assembly. This process is driven by complex and poorly understood atomic interactions of peptides with the lipid bilayer environment. The lack of suitable experimental techniques that can resolve these processes both at atomic resolution and nanosecond timescales has spurred the development of computational techniques. In this review, we summarize the significant progress achieved in the last few years in elucidating the partitioning of peptides into lipid bilayer membranes using atomic detail molecular dynamics simulations. Indeed, partitioning simulations can now provide a wealth of structural and dynamic information. Furthermore, we show that peptide-induced bilayer distortions, insertion pathways, transfer free energies, and kinetic insertion barriers are now accurate enough to complement experiments. Further advances in simulation methods and force field parameter accuracy promise to turn molecular dynamics simulations into a powerful tool for investigating a wide range of membrane active peptide phenomena

Ion dynamics in cationic lipid bilayer systems in saline solutions

Positively charged lipid bilayer systems are a promising class of nonviral vectors for safe and efficient gene and drug delivery. Detailed understanding of these systems is therefore not only of fundamental but also of practical biomedical interest. Here, we study bilayers comprising a binary mixture of cationic dimyristoyltrimethylammoniumpropane (DMTAP) and zwitterionic (neutral) dimyristoylphosphatidylcholine (DMPC) lipids. Using atomistic molecular dynamics simulations, we address the effects of bilayer composition (cationic to zwitterionic lipid fraction) and of NaCl electrolyte concentration on the dynamical properties of these cationic lipid bilayer systems. We find that, despite the fact that DMPCs form complexes via Na+ ions that bind to the lipid carbonyl oxygens, NaCl concentration has a rather minute effect on lipid diffusion. We also find the dynamics of Cl- and Na+ ions at the water-membrane interface to differ qualitatively. Cl- ions have well-defined characteristic residence times of nanosecond scale. In contrast, the binding of Na+ ions to the carbonyl region appears to lack a characteristic time scale, as the residence time distributions displayed power-law features. As to lateral dynamics, the diffusion of Na+ ions within the water-membrane interface consists of two qualitatively different modes of motion: very slow diffusion when ions are bound to DMPC, punctuated by fast rapid jumps when detached from the lipids. Overall, the prolonged dynamics of the Na+ ions are concluded to be interesting for the physics of the whole membrane, especially considering its interaction dynamics with charged macromolecular surfaces